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The study of hidden symmetries within Dirac’s formalism does not possess a 
systematic procedure due to the lack of first-class constraints to act as symmetry 
generators. On the other hand, in the Faddeev-Jackiw approach, gauge and 
reparametrization symmetries are generated by the null eigenvectors of the 
sympletic matrix and not by constraints, suggesting the possibility of dealing 
systematically with hidden symmetries throughout this formalism. It is shown 
in this paper that indeed hidden symmetries of noninvariant or gauge fixed 
systems are equally well described by null eigenvectors of the sympletic matrix, 
just as the explicit invariances. The Faddeev-Jackiw approach therefore provide 
a systematic algorithm for treating all sorts of symmetries in an unified way. 
This technique is illustrated here by the SL(2,R) affine Lie algebra of the 2- 
D induced gravity proposed by Polyakov, which is a hidden symmetry in the 
canonical approach of constrained systems via Dirac’s method, after conformal 
and reparametrization invariances have been fixed. 
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1 Introduction 


Invariant theories, both of gauge and reparametrization type, have been dealt with in 
the last 50 years by Dirac’s Hamiltonian formalism J]J, with its categorization of con¬ 
straints as first and second-class and the separation of equalities as weak and strong. 
Within this context the symmetries of a given theory are generated by a complete set 
of first-class constraints obtained by a well defined algorithm, and these constraints 
are used to reduce the original Hilbert space into the subspace of physical solutions. 
Nevertheless, noninvariant theories and gauge-fixed invariant models presenting only 
second-class constraints may still possess a set of (hidden) symmetries that cannot be 
unveiled by Dirac’s methodology. As we know, second-class constraints are used to al¬ 
ter the original (Poisson) bracket operation so that they can be imposed as operators 
identities, but cannot be used as symmetry generators. Within Dirac’s formalism, 
the disclosure of those symmetries that may be (possibly) hidden in a model will 
depend then crucially on our ability and experience on the subject, due to the lack 
of a systematic methodology. 

Recently, Faddeev and Jackiw|2|] proposed an alternative approach to constrained 
systems that avoids the separation of constraints into first and second-class and the 
use of weak and strong equalities. This new method of analysis has been successfully 
used by many authors|3], [I]. []. [| % || , and is by now a standard technique to deal with 
constrained systems. In a series of papers MM 0-0, this author and collaborators 
have proposed an algorithm, based on the Faddeev-Jackiw technique, to treat systems 
with constraints, both invariant and noninvariant, on the same foot: for noninvariant 
systems the sympletic matrix defining the geometrical structure of the theory will 
always end up being invertible even if it was not to begin with. This matrix inverse 
will then provide the Dirac brackets used in the quantization process(up to ordering 
ambiguities). For gauge theories, the sympletic matrix is never invertible, but it has 
been shown [0, |T2| that the null eigenvectors of this matrix are the generators of the 
intrinsic gauge symmetry of the model. 

It is our intention in the present work to show that those symmetries that persist 
after gauge fixing and are called hidden in the literature due to the lack of a set of first- 
class constraints to generate them, need no special treatment in the Faddeev-Jackiw 
context, being described by a set of null eigenvectors of the sympletic matrix exactly 
as for the explicit symmetries. While for those systems where Dirac’s formalism is 
free of pathologies, the advantage of Faddeev-Jackiw’s method over Dirac’s seems 
to be only of practical nature (it is simpler and quicker to use), the case of hidden 
symmetries shows a definite edge in favor of the Faddeev-Jackiw geometric approach. 
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To illustrate this systematic process of finding the hidden symmetries of an arbi¬ 
trary action, we study the SL(2,R) symmetry of Polyakov’s 2D induced gravity. The 
present status of this theory is as follows. The SL(2,R) current algebra symmetry of 
the induced gravity was discovered by Polyakov in a few years ago. Later on this 
symmetry has led to an exact solution of the 2-D quantum theory of gravity interact¬ 
ing with conformal matter fL|. It allows us to write down recursion relations for the 


(Euclidean) Green functions and to determine the renormalization of the parameters 
of the theory. This symmetry, which is hidden under the more conventional approach 
to current algebra, was discovered by Polyakov by analyzing the anomaly equation. 
In spite of its importance for the development of this very active field, the physical 
origin of this symmetry remained obscurcQ (for a nice review on the development of 
this theory and a complete list of references see Ref. ||1 5|| ). 


The study of this symmetry under the canonical point of view has been pursued 
by many authors 0 ^ 0. 0, mostly based on Dirac’s theory of constrained 
systems]!]]. The puzzling point, as indicated by Egorian and Manvelian|[T7|] and by 
Barcelos-Neto S. is that the local form of Polyakov’s model in the light cone gauge 
has only second-class constraints when written in Dirac’s light-front coordinates, and 
no constraints whatsoever if written in instant-front variables^]. Dirac’s method for 


dealing with constraints is well established and the basic feature of this approach is 
that gauge and reparametrization symmetries are related with first-class constraints, 
not present in this model. In [|17]], the (hidden) SL(2,R) symmetry was found following 
a somewhat modified Dirac’s method, but their solution was criticized in because 
they did not find the (supposedly also hidden) first-class constraints that generate this 
symmetry. This task has been pursued in |T9[ more recently. Working with the theory 
compactified in the (light-cone) space-dimension ( x ~) to a circle, the author was able 
to identify among the (infinite) set of second-class constraints, three apparently first- 
class constraints. He identified these constraints as the zero-modes of the set of 
second-class constraints of the gauge fixed theory. Unfortunately these constraints 
do not seem do correspond to the generators of the SL(2,R) symmetry, as mentioned 
by the author. This happens because these zero-mode constraints are in fact second- 
class, which explains why they fail to generate the required symmetry. 

An alternative approach to the problem has appeared in the literature quite re¬ 
cently |2T|. In this work, the authors were able to clearly identify the generators of 
this hidden symmetry as improper gauge constraints, as proposed by Regge and Teit- 
elboim [E2l |23], |24l] . The separation of first-class constraints into proper and improper 


^ee note added at the end of the paper. 
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was proposed in p2| in order to resolve some difficulties of Dirac’s method: some¬ 


times, among the second-class constraints there appears to exist a subset of first-class 
constraints. In P3 it was shown that a precise mathematical treatment of this subset 
of first-class constraints reveals that they are in fact second-class and, if treated as 
such, Dirac’s procedure may be consistently applied. Proper gauge transformations 
then, represent true symmetries of the theory and do not change the physical space of 
states. Improper gauge transformations, as defined in [E3, 23] do change the physical 


state of the system, mapping one physical solution into another different one. While 
proper gauge transformations can be eliminated by gauge fixing, improper gauge sym¬ 
metries cannot since that would exclude physically allowed solutions. Thus improper 
gauge transformations will remain as hidden symmetries. 

The outline of this paper is as follows. In section II we review the FJ formalism 
for invariant theories and show that the explicit symmetries, as well as the hidden, 
are generated by the null eigenvectors of the sympletic matrix, giving in this way an 
unified treatment for the subject under this approach. We end up this section with a 
discussion of two simple examples of physical systems presenting hidden symmetries. 
We show next, in Section III, that the sympletic matrix for the Polyakov model in 
the light cone gauge is singular and possess three null eigenvectors that happen to be 
the generators of the SL(2,R) residual symmetry which, therefore, becomes manifest 
in the Faddeev-Jackiw formalism. 


2 Faddeev-Jackiw Formalism For Invariant Theo¬ 
ries 

Faddeev and Jackiw|2|] showed that the treatment of a Lagrangian that is first-order 
in time derivatives dispense the use Dirac’s formalism. Consider a dynamical system 
with N bosonic degrees of freedom g* described by the following Lagrangian^] 

L — ai(q)qi — V(q) ; i = l,...,N (1) 

V(q) is the sympletic potential and dj(g) are the components of the canonical one- 
form a(q ) = di(q)dqi. This presentation of the theory is called non-standard since its 
equations of motion do not involve accelerations, and read 

5 We only consider a finite number of bosonic coordinates, the extension to the fermionic case ||2E|| 
and to field theory being straightforward]^]. 
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( 2 ) 



with the sympletic matrix f %3 being dehned by 


- ddj ddi 
= % “ % 


( 3 ) 


If this matrix is nonsingular, its inverse elements give the Dirac brackets of the theory 
([]]), with the dynamics being described by 


{F,G} 


dF 

dqi 



dG 

d qj 


( 4 ) 


However, when is singular the scheme above no longer works due to the existence 
of true constraints that reduce the number of independent degrees of freedom of 
the theory. These (sympletic) constraints will appear in this formalism as algebraic 
relations Q a (q) needed to maintain the consistency of the equations of motion upon 
multiplication by the (left) zero-modes vf of f t] . Then, from ([|) we get 


n a 


>i) T 


dV 

dqi 


where T stands for matrix transposition, and [] 


0 


( 5 ) 


Kfk = 0 ; a = 1,M (6) 

define the M independent (left) zero-modes of f tJ . 

There are basically three situations that can arise out of the scalar product defin¬ 
ing the sympletic constraints ([5]), which depends uniquely on the angle between v^ 
and (4- for each value of a. In the case of noninvariant systems all zero-modes are 
nonorthogonal to the gradient of the sympletic potential and constitute a true set of 
(sympletic) constraints that can be iteratively implemented into the canonical sec¬ 
tor of the original theory fll]), producing, step by step, a deformation of the original 
sympletic matrix, which end up being nonsingular ||10|| . In Dirac’s terminology, this 
corresponds to a theory with purely second-class constraints. The sympletic algo¬ 
rithm described above has been shown to lead to the same results as that of Dirac’s, 

6 For the case of field theory, where the zero-modes are in general operator valued objects, trans¬ 
position is defined in the sense of f [v T A)B = f A(vB). 
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although in a more economical way. A completely different situation appears when 
the gradient of the sympletic potential happens to be orthogonal to all zero-modes. 
Then condition (^J) vanishes identically, the equations of motion are automatically 
validated and no sympletic constraints appear. This happens due to symmetries in 
the sympletic potential, a situation typical of gauge and reparametrization invari¬ 
ant theories. There is no hope to deform the sympletic matrix into an invertible 
one unless a gauge-fixing term is introduced to break the symmetry Jll|, or one has 
the constraint solved and pass, via Darboux’s transformation, to a reduced canonical 
set[[2(J. Finally, we can have a mixed situation where some zero-modes are orthogonal 
to and some are not. In such cases one uses the sympletic algorithm to eliminate 
the constraints until only those zero-modes associated with symmetries remain. 

Let us assume, for definiteness, that all sympletic constraints have been elimi¬ 
nated, and therefore only the zero-modes associated with gauge symmetries are still 
present. In first-order such theory can be written as 


L = a k (q)q k + Va^ a {q) - V(q) 


( 7 ) 


where the sympletic variables q k \k = 1,..., N and r/ a ; a = 1,..., M form a set of 
gauge fields. We have also relabeled the Lagrange multipliers A a —> By 

hypothesis, the sympletic matrix (f km ) constructed out of the g fc -variables only is 
nonsingular, i.e., det / ^ 0. In terms of the q' s and rf s the sympletic matrix becomes, 


/ 


(/) 

Id fl\T 
y dq ) 



( 8 ) 


In this compact notation (/) represents a (NxN) nonsingular matrix defined in (|3|). 
(|h) represents a rectangular (NxM) matrix defined as 



ka 


( 

an 1 

an 2 

an M 

\ 

dqi 

dqi 

dqi 


an 1 

an 2 

an M 



dq 2 

dq 2 

dq 2 


\ 

an 1 

an 2 

an M 

) 

dq N 

dqN 

dq N 


( 9 ) 


The sympletic matrix (/) in (^) may have M zero-modes with the following (block) 
structure 


K) = 


( f 1 8Q. a 

\J km) 9qm 

l(“) 


( 10 ) 
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where the first element is an (Nxl) column and U“) is an (Mxl) column of zeros 
except in its a-th entry that is unity. For instance, 


l(2) 


( ° \ 
1 

0 


( 11 ) 


\°j 

Using (||) and O. the zero-mode condition ([5j) can be written as two sets of equa¬ 
tions; the first being automatically satisfied, while the second set is valid, as long as 
f2 a defines an algebra under /, i.e., 


( 12 ) 


(or ii*} = c^sr 

with the curly bracket operation being defined by (||). Therefore, the vectors (u(t) in 
(0) define a full set of null eigenvectors of the sympletic matrix (/) over the constraint 
surface. Since by hypothesis the gradient of the sympletic potential is orthogonal to 
all zero-modes, they must be the generators of the symmetry transformation that 
leave the action invariant. Specifically, we must have that the symmetry of the action 
over the constraint surface reads [pd] 


1 dQ a 

$eQm (/mn) ~ 


dq n 


Sjla = 


(13) 


The e a form a set of infinitesimal parameters that characterize the transformations. 

The point to be stressed here is that the symmetry transformations ( [hip may 
reflect either the gauge or reparametrization properties of an invariant theory, or 
they may as well be hidden symmetries of a noninvariant model. While the former 
is well described by Dirac’s formalism using the first-class constraints as symmetry 
generators, the later has not an easy description in the usual canonical approach. 

To give simple but nontrivial illustrations of how the sympletic formalism just 
described reveals a hidden symmetry, let us consider the examples of Floreanini- 
Jackiw chiral boson[p7]], and 2D Maxwell fields. The theory of self-dual bosonic fields 


proposed by Floreanini and Jackiw is described by the following action 


S= cPx(rf><fi — (ft'cp 1 ) 


( 14 ) 









where dot and prime have their usual meaning as time and space derivatives. The 
equations of motion possess a left mover held solution, as it should, but in addition 
to that there is a space independent freedom in terms of an arbitrary function of time 
h(t) as 


(j>(x, t) = cj) + (x + t) + h(t) (15) 

In terms of Dirac’s formalism, this theory has only a single second-class constraint, 
being therefore unable to disclose the origin of the symmetry (p~5|) . The only hint 
for its presence comes from the fact that Dirac’s constraint matrix does not have a 
unique inverse. On the other hand, the model’s sympletic matrix, 


f(x-y) = 2dj(x-y) (16) 

has as its only zero-mode, an arbitrary function of time which, according to (|T3|), 
accounts for the symmetry (ED above. To fix this invariance one has to adjust the 
boundary conditions such that the sympletic matrix zero eigenvalue equation has 
no nontrivial solutions which, in this case can be obtained with chiral boundary 
conditions. 

Another simple but nontrivial example is given by Maxwell electrodynamics in 
(1+1) space-time dimensions PS 


Following Dirac’s approach one is bound to find 


two first-class constraints 


Do — 7To 

0, = d x E (17) 

which are respectively the primary constraint and the secondary Gauss law constraint. 
Introducing the gauge fixing conditions 

Xo = A 

Xi = d x A (18) 

corresponding to the radiation gauge, the constraints above will form a second-class 
set. If the boundary conditions of the problem are chosen such that the Laplacian 
operator (d x ) has no nontrivial zero-eigenvalue solution, then the Dirac matrix of 
constraints will have a well defined and unique inverse. Certainly in such a gauge all 
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Dirac brackets vanish identically, indicating that no propagating modes survive after 
gauge fixing. This result could have been told in advance since the original phase- 
space having dimension four would result in a theory with zero degrees of freedom 
after four gauge conditions have been imposed, ffowever, if the boundary conditions 
are less restrictive, there will be some (hidden) residual symmetry left over which, as 
we will see, is responsible for the theta-vacua structure displayed by the model. Let 


us consider the physical system above enclosed in a finite space-time box[29 


S 



. _ TP TP^ 

4 ^ 


The sympletic matrix is easily computed to be 


(19) 


/O - y) 


l 0-1 0 \ 

1 0 -d x 5{x - y ) 

V 0 -d x 0 ) 


( 20 ) 


ffere the sympletic variables are £& = {A, E, A} with A — A±, A 0 = A and E = F 0 1 . 
From (H) and (|10|) we get the zero-mode of (|20|) as 


v(x) 



( 21 ) 


with u(x,t ) being a totally arbitrary function of the space-time coordinates. This 
zero-mode displays the well known Abelian gauge symmetry of the model which reads 


SA = d x u 

5E = 0 

SA 0 = d t u (22) 

After imposing the radiation gauge, the sympletic matrix assumes the form 


frad(% V ) 


( 0 

-1 

0 

d x 

\ 

1 

0 

-d x 

0 

0 

-d x 

0 

0 


\ d x 

0 

0 

0 

J 


5(x-y) 


This matrix displays the structure advanced in ([|) with 


(23) 
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and 


0 -1 

1 0 


S(x-y) 


(24) 


/ 0 


y) = 


'on' 


0 d x 
-d x 0 


S(x-y) 


and possess two zero-modes whose form is given by (IK);). For instance, 


(25) 



/ d x u \ 
0 
u 

V 0 ) 


(26) 


with a similar one for the second zero-mode. As discussed above, the existence of 
a sympletic matrix zero-mode in this case will depend on the Laplacian equation 
d 2 u = 0 possessing a nontrivial solution. If the boundary conditions are chosen such 
that this is the case, then we have that the two Laplacian eigenvectors 


Ml = Ci(t) 

u -2 = c 2 (t)x (27) 

do satisfy the following operator condition 

D^u k = 0 ; (no sum) (28) 

with 

DW = {d x ,xd x - 1} (29) 

It is easy to check that these operators close an algebra under the commutator oper¬ 
ation as 


[£>(*), = C* m D in) (30) 

where the only nonvanishing structure constants are C} 2 = C 21 = 1. One can see 
that redefining these operators as 
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( 31 ) 


D {1) £>« 

D {2) -> D {2) 


D {1) 

Dw 


the set D^ k> will form a simple representation of the Heisenberg-Weyl algebra. 

To see the physical consequences of this residual symmetry we follow Ref. p9 
closely. Recall that the solution for the Laplacian equation reads 


u(x,t) — a(t)x + b(t) (32) 

If one chooses u(t, j) = u(t, — j) = 0, then u(x,t ) vanishes identically and no sym¬ 
metry is obtained. However, if one chooses instead u(t, f) — u(t, — j) = 2n m, then 
u(t,x ) is nontrivial and is given by 


u(t,x ) = 2n^-N(t) 

R 

N(t ) = (n+ -n_)^ + i(n+ + n_) (33) 

with n± being defined at the time boundary by 


T 27 FT 

u ± (x) = u(x,±-) = —n ± 


(34) 


In [EUJ, the authors have called this solution at the time boundary as “classical vacua”. 
The zero-mode of the sympletic matrix therefore describes the interpolation between 
two topologically different vacua, u + (x) and U-(x), with relative winding number 
v = n + — ri-. Observe that the action for the Maxwell field now reduces to that of a 
particle moving on a circle 


S 



dt—N 2 

2 


l M (n + -n_) 2 
2 T 


(35) 


with the particle’s mass parameter being proportional to the circle’s inverse curvature 
M ~ if- Here 


v — 



(36) 


is the Pontryagin index and the zero-mode can safely be called as the instanton 
solution. 
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3 2D Induced Gravity 


In 2D there is no Einstein equation since the Hilbert-Einstein action is a topological 
invariant measuring the genus of the manifold on which one integrates. Nevertheless, 
conformally invariant matter fields induce an action over the metric holds through 
the conformal anomaly. Let us now consider the action for the induced 2D gravity 

m 


S = —k / y/^gRn 1 Rdrda 


(37) 


with k = c/967t, and c being the central charge of the matter held. This effective 
action is nonlocal in the metric helds but the nonlocality can be removed introducing 
an auxiliary scalar held 0(r, a) JT3], BUI BI 


S — — — J \f--g [g a ^d a (f) dp(j) + oiRcf^j drda (38) 

where ck is a constant parameter. The coupling of the scalar curvature R with the 
auxiliary scalar held </> breaks the conformal invariance of the nonlocal model. To 
break the reparametrization invariance we must impose some gauge conditions. We 
choose to work in the light cone gauge characterized by g /w = {g++ = 2h; g + _ = 
— 1; g _= 0}. With this choice the scalar curvature becomes 


R = -2 d 2 x -h 

Using these results in (|38|) we obtain, after two integrations by parts 


(39) 


L = dx 


'3 ,1 o V o + hio, o d x -<j> + <V) Of 


(40) 


As mentioned above, when described in instant-front variables[20], this model 


presents no constraints, while in light-cone variables only second-class constraints are 
present. This rehects the fact that all symmetries of the model have already been 
hxed. So, under Dirac’s approach, this model should present no symmetries. To 
examine the existence of a (hidden) symmetry under the Faddeev-Jackiw point of 
view we have to construct the sympletic matrix for this model. To this end we first 
relabel the Lagrange multiplier held as h —> d x +X, so that 


L = dx 


d x +(f) d x -(f) + d x +X(d x -<f> d x -<t> + oid 2 x -4>) 


(41) 


The sympletic matrix (in terms of (j) and A) reads 
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f(z,y) 


S(x-y) (42) 


( —2d x - adl--2dl-4>-2d x -(l)d x - 

\ -adl- - 2d x -4>d x - 0 

One immediately notices the similar structure of this matrix with that in (]^). There¬ 
fore, according to the discussion above, if there exist solutions to the matrix equation 

J dy~f(x,y)v {a \y) = 0 ; a = l,...,M (43) 

then they must corresponds to the zero-modes of f(x, y) that will reveal the symmetry 
remaining in the model. According to m the zero-modes must have the following 
structure 


I 5 dz e(y - z ) [ ad 2 z .g a {z ) - 2d z -(g a {z)d z -c 

9a(y) 


v^\y) = ^ 4 

which, after simple algebra, can be written in the following form 

v (a \y) = 


)( a ) (iA — ( \ ad y-9a{y) - g a {y)d y -(j) 
9a{y ) 


(44) 


(45) 


with g a (x) being a set of some yet undetermined functions. In components, the matrix 
equation (|45|) becomes the following pair of equations 


0 = J dy~d x -5(x~ - y~) (ad y -g a (y) - 2g a (y)d y -4>) 

+ J dx~ (ad 2 x - ~ ~ 2d x -(pd x -) 5(x~ - y~)g a (y) 

0 = f dy~ (—2d x -<i>d x - + ad%-) 5(x~ - y~) d y -g a (y ) - g a (y)d y -(pj 

It is a simple algebra to check that the first equation is trivially satisfied, therefore 
imposing no restriction over the functional form of g a {x). The second equation, on 
the other hand, is only satisfied if 

a 2 dl-g a (x ) = 0;a^0 (46) 

There are therefore three distinct null eigenvalues for the sympletic matrix (|42j) that 
reads 
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9 i(x) = e (1) (a: + ) 

g- 2 (x) = e < - 2 \x + )x~ 

g 3 (x) = e i3 \x + )(x~) 2 (47) 

If one removes these three zero-modes (with properly chosen boundary conditions) 

then the sympletic matrix becomes nonsingular and the elements of the inverse give 
the Dirac brackets that quantize the theory. Observe that the three differential op¬ 
erators satisfying 


D^g a = 0 ; (no sum) (48) 

are given by 

D {a) = {4-; x~d x - - 1; ( x~) 2 d x - - 2 x~j (49) 

that can be easily verified to be the generators of the SL(2,R) group algebra. Also, 
from ([L3|) we obtain the symmetry transformation of the scalar field under this group 
of transformations as 


S e (j) = — (V 1 ^" 1 ") + x e^ 2 \x + ) + {x ) 2 e ( ' 3 \x + )j d x -(j) 

+ ^ (e^(x + ) + 2oTe ( ' 3 ' ) (a; + )) (50) 

In conclusion, under the geometric approach of Faddeev-Jackiw, the SL(2,R) 
residual symmetry of Polyakov’s model manifest itself in terms of the zero-modes 
of the sympletic matrix, just as any other explicit symmetry, such as gauge and 
reparametrization invariances. While the study of symmetries inside Dirac’s formal¬ 
ism is based on the existence of first-class constraints, in the Faddeev-Jackiw formal¬ 
ism it relies on the orthogonality of the sympletic matrix zero-modes with the gradient 
of the (sympletic) potential. These (null) eigenvectors generate thus the isopotential 
lines of the theory, i.e. its symmetry lines, bypassing the first-class constraints for 
doing this job. This characterization of the Faddeev-Jackiw method in terms of null 
eigenvectors of the sympletic matrix, instead of constraints seems to be a definite 
advantage of this method over Dirac’s formalism since it avoids the use of first-class 
constraints, which in cases as the one above are certainly not available. 
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Note Added: After the completion of our work we became aware of the Ref.|32|, |33], |34| 
where these authors were able to recover Polyakov’s 2D gravity using the coadjoint 
orbit method. In particular Alekseev and Shatashvili |p2|| have shown the existence of 
a SL(2,R) symmetry in the 2D gravity by relaxing the periodicity boundary condition. 
The main point of this method is to exploit the fact that a coadjoint orbit of a group 
admits a sympletic two-form that may be used to construct an action that enjoys at 
least as much symmetry as the group whose coadjoint orbit is under consideration. It 
should be noticed that the method of our paper, using the Faddeev-Jackiw formalism, 
follows the opposite pathway, in the sense that we use the zero-modes of a singular 
sympletic two-form to determine the symmetries enjoyed by the action that produced 
that particular singular sympletic form. It is also worth to point out that the coadjoint 
orbit method is a highy mathematical subject while our method is acessible to all. 
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